How to make a fragile network robust and vice versa 
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We investigate topologically biased failure in scale-free networks with degree distribution P{k) oc 
fc""^. The probability p that an edge remains intact is assumed to depend on the degree k of adjacent 
nodes i and j through pij oc {kikj)~" . By varying the exponent a, we interpolate between random 
(a = 0) and systematic failure. For a > (< 0) the most (least) connected nodes are depreciated 
first. This topological bias introduces a characteristic scale in P{k) of the depreciated network, 
marking a crossover between two distinct power laws. The critical percolation threshold, at which 
global connectivity is lost, depends both on 7 and on a. As a consequence, network robustness or 
fragility can be controlled through fine tuning of the topological bias in the failure process. 

PACS numbers: 64.60.aq, 89.75.Hc 64.60.ah, 
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Scale-free networks, with power-law degree distribu- 
tion -P(fc) oc k~'^ , are remarkably resistant to random 
failure [J 0]. This quality is important when failure is 
to be avoided, as in the air-transportation network. It 
has been speculated that, also in nature, scale-free de- 
sign evolves as a way to achieve robustness % . On the 
other hand, robustness may be a problem when one tries 
to halt an epidemic. The fundamental question we ask 
and answer in this work is how one can delicately control 
whether a network is fragile or robust. 

Previous work has mostly concentrated on homoge- 
neous networks, in which all edges have the same chance 
to fail. However, by design or evolution the most critical 
edges of the network may become less prone to failure. 
Also, a targeted attack can disrupt the network after only 
a small fraction of edges fail [l|, . This shows that in 
heterogeneous networks the topology alone does not de- 
termine the susceptibility of breakdown. 

The critical properties of static phenomena and dy- 
namical processes are affected by the topology of the net- 
work of interactions ^5|]. It was recently shown [1,0] that, 
by accounting for a topology dependence in the interac- 
tion strength between the nodes, Jij oc (kikj)~°', one ob- 
tains a critical behavior that mimics the case of homoge- 
neous interaction but with a different degree distribution. 
The system with exponent 7 and topology-dependent in- 
teractions can be mapped to a homogeneous one, a — 0, 
but with an effective exponent 7', given by @, 0| 



/ 7 - a 
1 — a 



(1) 



We focus on failure in scale- free complex networks, me- 
diated by a dynamical process that depends on the net- 
work topology. Disregarding the presence of correlations, 
any such dependence has to be related only to the node 



degree fc. We have to choose between two possible ap- 
proaches, namely failure of the nodes or, as we implement 
here, failure of the edges. We express the failure probabil- 
ity for an edge between nodes i and j as qij = qij{ki, kj). 
We assume that the network depreciation occurs through 
a probability of occupation of an edge pij = 1 — qij , which 
depends on the degrees of the vertices. 



Pij oc Wij = {hkjY 



(2) 
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where Wij is the topology-dependent weight of the edge. 
Equation ([2]) is in the same spirit as the degree-dependent 
interaction proposed in [1, Since failure can be re- 
lated to the purely geometrical model of percolation, its 
understanding does not require interactions but can be 
achieved directly in terms of topological properties. 

A topology-dependent depreciation allows to interpo- 
late smoothly between random failure (a — 0) and inten- 
tional attack of links between hubs (a > 0), or intentional 
depreciation of edges between the least connected nodes 
{a < 0). We shall see that a 6 [2 — 7,1] defines the 
useful range of topological bias in the context of scale- 
free networks with finite mean degree (7 > 2). Degree- 
dependent failure was also studied in [8] for the case of 
node removal, while edge removal was investigated in . 

For uncorrelated networks, homogeneous random fail- 
ure (a = 0) can be solved using a mean-field ap- 
proach 0, [Hj [IHj [12] • Close to the critical fraction of 
occupied edges fc, the size of the largest connected clus- 
ter grows as (/ — /c)'', where the critical exponent (3 
depends on 7 fl^. For 2 < 7 < 3 the critical point 
vanishes, fc — > 0. As a consequence, we may say that 
networks with 7 < 3 are robust while networks with 
7 > 3 are fragile. All these results, however, are only 
relevant for the case of random failure. Henceforth, we 
will call the regime a > "centrally biased" (CB). The 
converse regime, a < 0, will be termed "peripherally bi- 
ased" (PB). 

To build our scale- free networks, we use the configu- 
ration model (CM) [ll|. The parameters of this model 
are the exponent 7, the number of nodes (or vertices) Ny, 
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FIG. 1: Average degree reduction in the depreciated network 
as a function of the original node degree. The main panel 
shows results for networks with 7 = 4 submitted to PB with 
a = ~1. One can see that the fraction k' /k grows as a power 
law with degree k, saturating at 1 at a scale that depends on 
the fraction / = n/Ne of edges in the depreciated network. 
The dotted lines are the numerical results obtained from 10 
network realizations of size iV„ = lO'^. For each one of these 
networks, the depreciation has been applied 10 times. The 
continuous lines are the best fit to the data of Eq. Q with 
the free parameter C(/) = 0.025, 0.22 and 0.92 for / = 0.1, 0.3 
and 0.9, respectively. The inset shows the same as in the 
main panel, but for a network with 7 = 2.5 subjected to CB 
with a = 0.5. In this case, the fraction fc'/fc decays with k. 
For small scales one may find a saturation depending on the 
fraction of edges in the depreciated network. The continuous 
lines are the best fits to the data of Eq. Q with values C(/) — 
0.28, 2.5 and 19 for / = 0.1, 0.3 and 0.9, respectively. 

and the minimum degree allowed kmin ■ Unless said other- 
wise, all the networks studied in this work have kmin = 2. 
Depending on these parameters and on the particular re- 
alization, we obtain a different number of edges N^. In 
this model, the degrees of the nodes are determined ini- 
tially from the desired distribution and then connections 
are assigned at random. 

To study the depreciation process, we start from a to- 
tal failure scenario, i.e., with all edges being initially re- 
moved from the network. We then gradually include the 
edges back, with probability proportional to some weight 
Wij that we will assume to follow Eq. ([2]). By stopping 
this process at intermediate steps, we can obtain results 
for the percolation problem as the fraction of occupied 
edges grows from zero to one. 

We now determine the probability Pij{f) that a par- 
ticular edge connecting nodes i and j is present in the 
network after a fraction / of the edges has already been 
included. This probability can be identified as Pij{f) — 
1 — Y\t=i ~ "^ij /■^t)i where Zt is the mean (over the 
inclusion process) of the sum of weights of all edges that 
have not yet been included in the network at step t, and 
n = fNf, is the number of included edges. Assuming 
Wij <^ Zt, we can write, 

P,,(/) «l-e-^(^')('=''=^)"°, (3) 
where the parameter D{f) can be determined using 



Y^Pij=n. Using wij as defined in Eq. ([2]), the 
Kasteleyn-Fortuin construction [l3| allows us to draw 
a parallel between the probability Pij and the degree- 
dependent interaction previously proposed in [6^, ^ . 

We now define Pk{f) as the mean probability that an 
edge from a node with degree k is present in the de- 
preciated network. We then have to average Pij over 
the nearest neighbors of a node with degree k, pk = 
C„„ ^"(^") (1 - \~D{f){kkn)-'^]) dfc„. For uncor- 
related networks the degree distribution of a neighbor is 
given by P„(A:„) = P(fc„)fc„/(A:). Performing the inte- 
gration and examining the asymptotic behavior of the 
resulting incomplete Gamma function, wc find that pk is 
well approximated by 

Pfe«l-e-^(«'=-°, (4) 

with C{f) = t^^^Dif), provided a G [2-7,1]. 
The same range of a is also featured in previous work 
on networks with degree-dependent interactions [1, 0]- 
Equation ([4]) is confirmed by the numerical results shown 
in Fig. m 

Equation. ^ can be used to determine the average 
degree of a node after depreciation, k'(k) — k pk- From 
that we can obtain the degree distribution of the de- 
preciated network: P'{k') = k'^ for C{f)k-°' > 1, 

and P'{k') = k''^ for C{f)k^" -C 1, where 7' is given 
by Eq. ([1]). We find that the degree distribution af- 
ter depreciation exhibits a crossover at a scale given by 
kg oc C{f y^". As expected, the crossover is not present 
in the random failure case, a = 0. However, if the failure 
process is affected by the topological properties of the 
network, as modeled by Eq. we have a characteristic 
scale kg that has not been observed before. The presence 
of this crossover is supported by the numerical results 
shown in Fig. O It is interesting to note that whether 
7 or 7' controls the decay at large degree depends on 
the sign of a. If a > (CB), we have 7 < 7' and 7' 
controls the decay at large degree, while for a < (PB) 
the larger exponent, 7, is the controlling one. Thus, the 
largest of the two exponents 7 and 7' controls the asymp- 
totic decay. A robust network with 7 < 3 under CB and 
a fragile network with 7 > 3 under PB may result in 
networks with similar degree distributions after depreci- 
ation. The numerical results shown in Fig. [5] correspond 
to the degree distributions of networks under CB and PB 
failure. 

Next we investigate the critical behavior associated 
with percolation. A network is above the critical point 
when a node connected to another node in the span- 
ning cluster has on average at least one other connec- 
tion, thus assuring that the cluster does not fragment. 
For an uncorrelated network, this condition is equivalent 
to (fc'')/(fc') > 2 0. In order to determine the critical 
fraction fc we perform the depreciation process until this 
critical condition is reached. In the simplest case where 
a = 0, if 7 > 3, the ratio {k''^)/{k') converges to a finite 
value as A'^i, — > cxo. In this case /c > 0, characterizing a 
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FIG. 2; Degree distribution of the depreciated network. 
This result was obtained for 10 network realizations of size 
Nv = 10^. For each network realization, the depreciation 
has been applied 10 times. In the main panel we show the 
degree distribution for networks with 7 — 2.5 submitted to 
CB failure with a = 0.5. From Eq. ([TJ we expect the value 
7' = 4.0 for the depreciated network. As a guide to the eye, 
the dotted lines indicate the power-law decays with exponents 
2.5 and 4. One can see that for small degree the distribution 
initially decays with a slope very close to 7, and then crosses 
over to a decay with a slope close to 7' at a scale that de- 
pends on the fraction / of edges in the depreciated network. 
This shows that the topology-dependent failure process intro- 
duces a characteristic scale in the degree distribution of the 
originally scale-free network. In the inset we show results for 
networks with 7 = 4 subjected to PB failure with a — ~1, 
which is equivalent to an exponent 7' = 2.5. Contrary to the 
result of the main panel, one sees a crossover from a slope 
7' at small degree to a slope 7 at large degree. The mini- 
mum degree was set to kmin = 4 in the results of the inset 
and kmin = 2 in the main panel. Surprisingly, for the case 
/ = 0.5 the degree distributions in the inset and the main 
panel are remarkably similar. To illustrate this similarity we 
included the results for / = 0.5 from the main panel in the 
inset and vice versa; these are the dashed (blue) lines. Al- 
though the two cases start with distinct degree distributions 
at / = 1, we obtain similar distributions at a certain point of 
the depreciation process. 

fragile network regime. On the other hand, when 7 < 3, 
one obtains (fc'^) oc k"^^, where k^ax is the largest de- 
gree of a finite network. If no other constraint is imposed, 
k^a. OC resulting in oc N-^'-^^^^-'~'\ 

for 2 < 7 < 3. As long as 7 < 3, /c — > as A'^i, ^ 00, 
characterizing the robust network regime. 

For CB {a > 0) it is possible that 7 < 3 while 7' > 3. 
Since the tail of the distribution at large values of k' 

decays as P'{k') oc k'^^ , the second moment (/c'^) no 
longer diverges and a robust network becomes fragile un- 
der CB. On the other hand, for PB {a < 0), 7' < 7, 
and the larger exponent, 7, should control the decay of 
the tail of the degree distribution. Therefore, one may 
think that for 7 > 3, {k''^) also does not diverge under 
PB and that a fragile network cannot turn robust. This 
simple reasoning is mistaken, however, as we show in the 
following. 



It is possible that the crossover scale kg becomes so 
large that in practice it does not influence a finite net- 
work. That is the case, for instance, of the distribution 
for / = 0.1 and = 10^, shown in the inset of Fig. 
It may be the case that, at the critical point of PB fail- 
ure, a network with 7' < 3 never displays an observable 
crossover, irrespective of the system size, that is, 

hifc{N,)) > kmaxiN,)- (5) 

When Eq. (O holds true, Eq. ^ may be rewritten as k' = 
C{fc)k^ °'- In this limit one can find a linear relation 
between the parameter C{fc) and the occupation fraction 
= {k')/{k) = C{k^-°')/{k). The second moment can 
be identified as {k'^) = C^{{k^-^°') - {k^-^°')) + C{k^-°'). 
From the critical condition {k'"^) / {k') > 2 we get 

As long as 2 < 7' and 2 < 7, the moments (fc^~") and 
(fc) should both converge to finite values independent of 
N^. The moment {k^~'^°') may or may not diverge, but at 
large scales it will grow slower than (fc^^^Q^ fc^^a-T^ 

Thus, considering that kmax oc Ny^^'' we arrive at 

the behavior 

/e«iV7^. (7) 

This result shows not only that PB can turn a fragile 
network robust but also that the critical exponent with 
which the threshold fc approaches zero is the same as 
expected for normal percolation (a = 0) for a network 
with a degree distribution decaying as P{k) oc k~'^ . 

We can now check the self-consistency of our initial 
assumption, Eq. ([5]), that networks with 7' < 3 at the 
critical point of PB failure do not present a crossover. 
As mentioned before, the crossover scale is given by ks oc 
At the critical point Eq. ([7]) then implies kg oc 
^h'~3)/[aiY-i)]^ while k^ax OC N^^-^'^l From this we 
obtain (h/kmaxy"'^''-^^ oc iV^-f-". As long as a < 
3 — 7, the crossover scale grows faster than the maximum 
degree, implying that critical networks with 7' < 3 and 
sufficiently strong PB do not display a crossover in their 
degree distributions. However, for weak PB, 3 — 7 < a < 
0, Eq. ([5]) is violated and the network may remain fragile. 

Figure [3] shows numerical results confirming that a ro- 
bust network with 7 = 2.5 submitted to CB failure with 
a — 0.5 turns fragile. In contrast, the second set of re- 
sults demonstrates that a fragile network turns robust 
even for a = 3 — 7 = —1, which is on the borderline 
between weak and strong PB. 

Our assumption that the probability of failure depends 
on degree k can be justified in different contexts. In ar- 
tificial networks, e.g., air transportation [isj . the capac- 
ity of the nodes scales with fc. Depending on whether 
k or capacity grows faster, this system should be bet- 
ter modelled by CP or PB failure, respectively. Soft- 
ware systems and metabolic networks [iTj consist of 
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FIG. 3: The values /c of the fraction of edges in the de- 
preciated network at the critical condition as a function of 
the network size A^„. Criticality is defined as the point where 
{k'^)/{k') = 2 [2], To compute this critical fraction we average 
over 10* network realizations for each set of parameters. For 
each network we apply the percolation processes 100 times. 
For networks with 7 = 2.5 submitted to CB with an effec- 
tive value 7' = 4 (continuous black line), we observe that the 
critical fraction /c converges to a finite value as A'^i, grows, 
confirming the conjecture that a robust network may turn 
fragile under CB. The opposite case, a network with 7 = 4 
submitted to PB with an effective 7' = 2.5 (dashed red line), 
has a critical fraction that decays with A'^„ as a power law, 
fc oc Nv^^" . The best fit to the data in this case results 
in 1/i^ = 0.35 ± 0.02, consistent with the value 1/3 expected 
from Eq. (0. This result shows that a fragile network under 
PB can behave in the same fashion as a robust network with 
a degree distribution controlled by 7' under random failure. 



many agents/nodes acting together in some function. If 
all agents are needed, the lack of any of them can inter- 
rupt the process. Alternatively, if any of the agents can 
start it, only removal of all edges halts the process. In 
both cases, depending on k, the edges turn more fragile 



or robust. Further, if disrupting the network is desirable, 
as in gene fusion networks of cancer development 18] or 
terrorist networks, the design of a dynamical process that 
targets links between the most connected nodes (CB) 
would be more efficient to globally break down the sys- 
tem. Also, to reduce the risk of epidemic spreading, it is 
better to disinfect/immunize connections between hubs 
than connections between small (air-)ports. Note that 
our analysis does not account for dynamical correlations 
in the failure process. It can happen that removal of a 
single edge triggers a breakdown, even if this edge only 
links to one of the least connected nodes [l^ . 

We conclude that topologically biased failure can have 
a dramatic effect on the percolation properties of scale- 
free networks. For central bias (CB, < a < 1), the 
degree distribution initially decays with the exponent 7 
up to a certain scale that depends on the fraction of oc- 
cupied edges, and then crosses over to a decay with an 
exponent 7' > 7 defined as in Eq. ([T]) . For peripheral bias 
(PB, 2 — 7 < a < 0) the crossover is also present but with 
7' controlling the early decay and the exponent 7 > 7' 
appearing at large degree. Our results also demonstrate 
that a robust network, for which the critical fraction /c 
converges to zero as the network grows, may turn frag- 
ile when subjected to CB [a > 0). Conversely, a fragile 
network, for which the critical point is larger than zero 
at any system size, may become robust when subjected 
to strong PB, a < 8 — 7. Fragility or robustness of a 
network is thus not only dependent on the exponent 7 
but can be tuned quantitatively by the exponent a char- 
acterizing the topological bias. 
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